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Abstract 

The synthetic aperture radar (SAR) technology enables satellites to efficiently ac- 
quire high quality images of the Earth surface. This generates a significant traffic from 
the satellite to the ground stations, and, thus, image downlinking may become a bot- 
tleneck in the efficiency of the whole system. In this paper we address the downlink 
scheduling problem for the Canada's Earth observing SAR satellite, RADARSAT-2. 
Being an applied problem, downlink scheduling is characterized with a number of con- 
straints that make it difficult not only to optimize the schedule but even to produce 
a feasible solution. We propose a schedule generation procedure that lets us nicely 
incorporate all the constraints and then effectively optimize the schedule. Our compu- 
tational experiments conducted on the real data show that the proposed algorithm is 
a significant improvement over the scheduling procedure currently in use. 
Keywords: Satellite, Scheduling, Optimization, Heuristics. 



1 Introduction 

Efficient scheduling of image acquisition and image downlinking plays a vital role in satellite 
mission planning. These operations are often interlinked and solved using scheduling heuris- 
tics that take advantage of the flexibility allowed within such integrated systems. Most of the 
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literature on satellite mission planning (image acquisition and downlinking) is concentrated 
on satellites with optical devices for image acquisition f2T\ or on modern satellites that 
use Synthetic Aperture Radar (SAR) technology jH El El |P| . We study the mission planning 
operations of Canada's Earth observing SAR satellite, RADARSAT-2. 

Our problem deals with scheduling image downlinks optimally where an image acquisition 
schedule is given. The acquisition of RADARSAT-2 images is scheduled by a separate 
algorithm based on client requests using information about available time windows, and 
it is not considered in our study. The downlink scheduling process currently in operation 
for the RADARSAT-2 mission makes use of a greedy type algorithm followed by human 
intervention where necessary. The operators of RADARSAT-2 are facing increased demand 
for its imagery and are exploring various methods to improve the efficiency of the downlink 
scheduling algorithms. We developed efficient heuristics to handle the downlink scheduling 
which resulted in significant improvements. Experimental analysis using historical data 
established that for situations where downlink scheduling is oversubscribed a reduction in 
unsuccessful downlink requests of 20 to 30 percent is possible using our algorithm. It should 
be noted that the unsuccessful downlinks observed in the test data set are almost all for 
image requests that are considered low priority 'background' acquisitions. For test data 
sets from other seasons with lower order volumes our algorithm produced schedules with 
zero unsuccessful downlink requests whereas the existing model experienced a loss rate of 
5 to 10 percent of low priority acquisitions. In all seasons, the total latency of downlinks 
is also reduced. These impressive experimental results were achieved using very large scale 
neighborhood (VLSN) search techniques [H [2] along with special ejection chain schemes [H] 
to design the associated neighborhood structure. 

The satellite image downlink scheduling problem (SIDSP) and its variations have been 
studied by many authors. Some of these works were concentrated on case studies for specific 
satellites [9] whereas other are more general purpose in nature [8| ITO t [TT | [T2 | [T3 | IT8 | l20 t [2T | 122] . 
Literature from machine scheduling [6l [17] and resource-constrained project scheduling [15] 
are also relevant in solving the SIDSP. However, each mission planning problem has its own 
restrictions and properties that can be exploited. Special care is needed to make sure that 
such restrictions are handled adequately which sometimes changes the inherent combinatorial 
structure of the problem significantly. Thus, investigating case studies of special SIDSPs are 
interesting and relevant as established in this study, although existing literature on the 
SIDSP considerably infiuenced our work. 

The paper is organized as follows. In Section m we give a brief description of our SDISP 
and introduce various notations and definitions. Section [3] deals with our heuristic algo- 
rithms. Data analysis and generation of test instances are reported in Section H] followed by 
computational results in Section |5] and concluding remarks in Section |6l 

2 Problem Description 

The SIDSP is defined by one satellite that orbits the Earth to acquire images and a set G 
of stationary ground stations that receive images from the satellite for further processing. 
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The image acquisition schedule is generated beforehand. Let Rhe a set of n image downhnk 
requests that needs to be downUnked to a ground station during the planning horizon of 24 
hours. For each request r E R, a release time br, deadline e^., downlink duration dr, priority 
Pr and a ground station Qr & G are prescribed. The interval [6^, e,.] is called the time window 
of request r. 

It may be noted that image downlinking needs to satisfy a variety of constraints and it 
may not be possible to downlink all the available images and, thus, rejection of requests is 
allowed. An image that is not downlinked before its deadline is considered ^ unscheduled\ 
The SIDSP deals with finding an image downlink schedule so that an appropriate utility 
function is maximized. The utility function considers the number of downlinks scheduled, 
their priority values, and the earliness of the downlink start time relative to its time window. 
The SIDSP is NP-hard since several NP-hard machine scheduling problems arc special cases 
of it. Hereafter, the terminology satellite refers to the spacecraft RADARSAT-2, which is 
the focus of this study. 

A downlink activity can be carried out only when the satellite is passing over the ground 
station coverage area, called station visibility mask, of the station. The satellite has two 
antennas that are used for downlinking and each ground station g E G has one or two 
channels that are used for receiving the downlinked image. If the ground station has two 
channels, two downlinks can happen in parallel (simultaneously) under certain circumstances 
which will be discussed later. In our case study, around 2/3 of the ground stations have one- 
channel capability and 1/3 of the stations have two-channel capability. Ground station 
channels work independently. Satellite antennas also work independently and are capable of 
downlinking to two different ground stations simultaneously. 

Ground stations are also classified based on their transmission power. Let Gi be the 
set of ground stations that are in half-power setting and G2 be the set of ground stations 
that are in full power setting. Then G — G^U G2 and Gif] G2 — $■ When two images are 
downlinked one after another to stations with the same power setting, there must be a gap 
of S seconds between the two downlinks. When two images are downlinked consecutively to 
two stations with different power settings, the required gap between the downlinks is A > 5 
seconds. For full power ground stations, only one downlink is allowed at a time. 

There are two ways of image downlinking: pass-through mode and store-forward mode. 
In pass-through mode, an image is downlinked as it is being acquired. This is possible only 
when the area to be imaged is within the visibility mask of the ground station to which it is 
downloaded. In store-forward mode, images taken earlier are stored in the satellite recorder 
are downlinked as opportunity arises. 

The visibihty mask of a ground station g E G can be represented by a collection Vg of 
non-overlapping time intervals. We call Vg the normal visibility mask. Certain downlink 
requests require better reliability and in such cases we use finer components to define a 
visibility mask Vg which we call high reliability visibility mask. Thus, for each g E G, we 
have normal and high reliability visibility masks. 

Let tr be the scheduled downlink start time of request r E R. li the request r is unsched- 
uled, the value of tr is undefined and this is denoted by tr — 0. The set S — {tr : r E R} 
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determines the solution to the SIDSP and it is called a schedule. Let R{S) — {r : re 
R and tr 7^ 0} be the set of scheduled requests under S. Recall that the priority of a 
request r E R is Pr- Thus, a natural objective to consider is to maximize '^reR{s)Pr- 

This objective measures the total priority level achieved by a downlink schedule. However, 
for the heuristics we considered, this objective produced poor solutions probably because 
the number of local maxima turned out to be very large. Since we are also interested 
in scheduling a request as early as possible, we have modified the objective function to 

be f{S) — '^reR{S)Pr ^1 — • e ^Id^-b ) ' "^^cre < a < 1 is a parameter reflecting the 
importance of scheduling each request downlink as early as possible. Note that when a = 0, 
the objective function reduces to maximizing ^reR{S)Pr- Thus, our model for the SIDSP 
reduces to 




subject to 5* e F , 



where F is the collection of all schedules S satisfying the following constraints and assump- 
tions: 

(1) A downlink cannot start earlier than the release time of the corresponding request and 
must be finished by its deadline. If a downlink cannot be done within this time window, 
the request is unscheduled. 

(2) If there are more than one request to downlink one image to different ground stations, 
either all or none must be scheduled. 

(3) There must be a gap of at least 6 seconds between two consecutive downlinks under the 
same power setting. 

(4) The satellite antennas work either in half-power setting or full-power setting. The tran- 
sition from one power setting to another requires A > ^ seconds. During the transition, 
none of the two antennas can work. Thus, if two consecutive requests are scheduled to 
two different ground stations with different power settings, then there has to be a time 
gap of at least A seconds between these downhnks. 

(5) If a satellite antenna is in the full-power setting, then only this one antenna can work and 
the other has to be idle. For both of the sateUite antennas to be working independently, 
the antennas have to be in the half-power setting. An antenna in full power setting is 
used if and only if the corresponding ground station is in the full power setting. 

(6) We assume that the satellite is in the half-power setting at the beginning of the planning 
horizon and guarantee that it is in the half-power setting at the end. We also assume that 
downlinking (in half-power setting) can start right from the beginning of the planning 
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horizon, i.e., there was no downhnk activity for at least 5 seconds before the beginning 
of the planning horizon. 

(7) No downlink activity happens outside the planning horizon. 

(8) Once a downlink starts, it cannot be preempted, i.e., an image cannot be split. 

(9) A pass-through request is considered as urgent and, if possible, has to be scheduled at 
the earliest possible time in its time window. 

Note that our objective function does not explicitly distinguish between urgent and non- 
urgent requests. However, this is taken care off by running our algorithm in two phases, 
where in the first phase all urgent requests are scheduled followed by scheduling of non- 
urgent requests in the second phase. 

3 Development of the Algorithm 

Our solution approach is based on a local search with an ejection chain neighborhood. A 
high-level description of the algorithm is given as follows: 

Pre-processing: This step modifies the data to handle constraint ([3]). 

Urgent requests processing: This step schedules all the urgent requests using the ejec- 
tion chain algorithm. 

Non-urgent requests processing: Starting with the partial schedules of urgent requests, 
we schedule non-urgent requests using the ejection chain algorithm, while maintaining 
the scheduled urgent requests intact. 

3.1 Pre-processing 

We first perform a pre-processing step which modifies the data to simplify the presentation 
of the algorithm and efficiently handle some of the constraints. One of the constraints of 
the SIDSP is to have a gap of at least 6 units between two consecutive downlinks. To 
accommodate this, we add 6 to dr and for all requests r & R. Also, we increase by 6 the 
upper bounds of each of the time intervals in normal and high reliability visibility masks for 
all (/ e G. Thus, even if two downlinks are scheduled such that one is followed immediately 
by another, a gap of size S between the downlinks is guaranteed. 

Note that such an adjustment in the data prohibits some downlink activity at the end 
of the planning horizon which otherwise could have been scheduled. One can increase the 
end of the planning horizon by 6 to take care of this. In that case, the conflicts between 
downlinks at the end of the previous planning horizon and the beginning of the current 
planning horizon can be taken care off by using appropriate perturbations and careful data 
handling. 
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Further, we replace A by A — 5 because the gap between downhnks to the stations when 
different power settings are used will always include the additional 6 units of time we added. 
Hereafter, we assume that problem data is modified to reflect the changes discussed above 
and, hence, we can ignore explicit treatment of constraints ([3]). 

3.2 Schedule Generator 

This section describes the greedy algorithm we used to generate a schedule form a given 
ordered subset of requests R. Let us now discuss some terminology and notations to simplify 
the presentation of the algorithm and its validity proof. 

Let y be a finite set where elements of V are non-intersecting intervals. We call sets 
of the type V interval set. Since the elements of V are non-intersecting intervals, V can 
also be viewed as an ordered set with the natural order produced by the position of these 
intervals on the real line. Thus, V is represented as y = {[^i, U2], [^2, ^2], • • • , [^v, u^]} where 
ii < ui < £2 < U2 < ■ ■ ■ < iv < Uv and v = \V\. We refer to the zth interval in V as Vi. 

Consider an arbitrary interval [a,b]. Let us introduce sets of intervals J, L and R as 
follows: 

B = {[ij, Uj] eV : ij <b and Uj > a} , 

C = {[ij, a] : 3j such that [ij,Uj] G V and ij < a < Uj} and 

D = {[b, Uj] : 3j such that [ij,Uj] G V and ij < b < Uj} . 

Then subtracting the interval [a, b] from V generates the set V Q [a, b] of non-intersecting 
intervals given by 

VQ[a,b] = {V\B)UCUD. 

The notation V Q [a, b] is read as V minus interval [a, 6]. 

An interval [a, b] said to be a subinterval of the intervals set V if there exists k such that 
ik ^ CL < b < Uk and [ik, Uk] G V. This relationship is denoted by [a, b] d V. 

Let R* be an ordered subset of R. Given R*, we now present a greedy algorithm to sched- 
ule requests in R* following the order prescribed for R*. The algorithm maintains several 
indicator interval sets representing channel availability at ground stations and antenna avail- 
ability on the satellite. After scheduling a request, the algorithm updates these indicator 
sets. 

Let A be an interval set that represents the time intervals when both satellite antennas are 
available. Let H and F be interval sets indicating the time intervals when a half-power and 
full-power downlinks can happen, respectively. Let Qg and Qg be interval sets indicating the 
availability of the ground station G G in normal and high reliability visibility, respectively. 
Note that we need only one pair of indicator sets for a two channel ground station. Indeed, 
the number of channels in this case is not limiting since the number of simultaneous downlinks 
is constrained by the number of antennas. Hence, only the visibility of the ground station 
has to be considered. 

The greedy algorithm works as follows: 
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1. 



Let R** ^ R*. 



2. Initialize the sets A, H and F to {[0, 24 hours]} which is the planning horizon. Initialize 
the set Qg to the normal visibility mask of g and to the high reliability visibility 
mask of g for every ground station g E G. Note that for every lij] G there is an 
interval [Ik, Uk] G Qg such that £fc < < < Wfc. 

3. Select the first request r G i?**. 

4. Suppose that r is downlinked to a half power station g e Gi. Let Q ^ Qc, if requires 
normal reliability and Q otherwise. Choose the smallest i & {1,2, . . . ,\H\} and 
j G {1,2,...,|(5|} such that | ifj Pi Q j H [6^ , 6^]] > dr- If no such i and j exist, the request 
r is unscheduled and tr ^ 0. Otherwise tr x, where Hi fl Qj fl [6,., e^] = [x, y]. 

If tf. 7^ 0, the indicator sets are updated as follows. Let X = [tr,tr + dr]. Apply 
X XQ[ii, Ui] for every [ii, Ui] G A. Now X contains all the intervals inside [U, U + dr] 
during which exactly one antenna was available. Apply H ^ H Q [ii,Ui] for every 
[£i, Mi] G X. Also apply A ^ AQ [tr, U + dr] and F ^ F Q [U - A,tr + dr + A]. 
Finally, if is a one channel ground station, update Qg Qg Q [tr,tr + dr] and 

Q]^ QlQ[tr,tr + dr]. 

5. Suppose that r is downlinked to a full power station g E G2- Let Q ^ if r requires 
normal reliability and Q Q\ otherwise. Choose the smallest iG{l,2,...,|F|} and 
i G {1,2,...,|(5|} such that \FiP\Qj fl [6^, ^r] | > (i^. If no such i and j exist, the request 
r is unscheduled and tr ^ 0. Otherwise, tr ^ x, where Fj fl Qj fl [br, e^] = [x, y]. 

li tr 7^ 0, the indicator sets are updated as follows: F -(^ F Q [tr,tr + dr] and H <(— 
-ff [tr — ^7 + c^r + ^]- Note that it is not necessary to update the indicator set A 
since no downlink can happen if neither H nor F is available. 

6. Update R** to R** \ {r} and, if \R**\ > 0, proceed to StepO 

7. For every pair (r, u) of dual requests, if either tr = and 7^ or tr 7^ and = 0, 
remove both r and m from R* and proceed to Step El 

8. Assign the downlinks to the particular ground station channels and satellite antennas. 

Theorem 1. The greedy algorithm terminates with a feasible schedule of requests in R* in 
0{n^) time. 

Proof. Our updating scheme of the intervals sets A, H and F ensures that no antenna 
conflict arises and all the downlinks happen within the planning horizon. By subtracting 
\fr — ^,tr + (ir + A] from F for every half power downlink request r, we guarantee that 
no full power downlink will happen within A units from the downlink r. Similarly, no half 
power downlink can happen within A units from from a full power downlink r. Also, the 
updating of the indicator sets Qg and guarantees that no channel conflicts occurs and the 
downlinks obey visibility mask constraints. The pre-processing of data assures that there is 
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a gap of at least S units between two consecutive downlinks. Finally, dual requests constraint 
is satisfied due to Step [7] of the algorithm. Thus, the schedule generated is feasible. 

Let us now analyse the complexity of the algorithm. The primary operations in each 
iteration is: (1) to find the smallest i and j to satisfy certain condition and (2) update the 
indicator sets A, H, F, Qg and Q^. Note that the size of the indicator sets Qg and Q^g for 
some g E G may increase but it is bounded by 0{ng) where Ug is the number of requests to 
be scheduled to the station g. Similarly, the sizes of the indicator sets A, H and F are limited 
by 0{n). Thus, operation (1) for these sets can be performed in 0{n) time by simultaneous 
scanning of the sets. Operation (2) can also be performed in 0{n) time. Indeed, we only 
need to update a fixed number of indicator sets, and updating each of them takes 0{n) time 
(for a downlink to a half-power setting ground station, manipulating with X also needs only 
0{n) time; note that each of the H ^ H Q [ii, Ui] for [ii, Ui] G X operations needs only 0(1) 
since all of these operations affect only one interval in H). The number of iterations is at 
most O(n^) and, thus, the algorithm complexity is O(n^). □ 

3.3 Construction Heuristics 

The greedy scheme discussed above can be used to obtain various construction heuristics 
for the SIDSP. Choose R* = R and depending on how we order elements of R* we obtain 
different construction heuristics. We considered sorting the requests based on one or more 
of the following fields: 

1. Priority pr] 

2. Time window range w* = [e.^ — br — dr\. Note that if we do not use rounding, time 
windows for all the requests could be different and, thus, sorting using a secondary 
criterion became somewhat superfluous. 

3. Downlink duration dr- 

After extensive computational experiments using various combinations of the above fields 
to sort requests, it is determined that sorting elements of R* with respect to pr first, ties 
are broken by sorting with respect to w* values and further ties, if any, are broken by dr 
values provides the best performance. The resulting solution is used as the initial solution 
for our improvement heuristic that uses ejection chain neighborhoods. This Ejection Chain 
algorithm is described in the next section. 

3.4 Ejection Chain Improvement Heuristic 

We now show that we can find improved solutions using ejection chains. Note that the quality 
of the solution produced by the greedy algorithm depends on the ordering of elements in 
R*. By repeating the algorithm with different orderings we get possibly different solutions 
and the overall best solution can be chosen. However, exhaustive searching over all such 
orderings is impractical. Our ejection-chain algorithm chooses 'good' ordering for R* so that 
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the output generated is of high quahty. Let p be the set of all permutations of elements of 
R. For any vr G let S.^ be the schedule produced by the greedy algorithm by choosing 
R* = R and the ordering of elements of R is given by n. Thus, we concentrate on solutions 
{Sn : TT G to design our ejection algorithm. In this sense, we consider p as the family of 
feasible solutions of the SIDSP. Thus, to solve the SIDSP, we essentially need to solve the 
downlink requests permutation problem (DRPP) given by: 

Maximize 0(7r) = /(S'^r) 
subject to TT G p. 

Ejection chain based neighborhoods are commonly used in developing VLSN search algo- 
rithms [H [2] for complex combinatorial optimization problems. For example, the well known 
Lin-Kernighan heuristic and its variations are ejection chain algorithms that are extremely ef- 
ficient for the travelling salesman problem [3]. We now use the idea of the ejection chains 
to develop a simple and effective heuristic to solve DRPP (and, hence, the SIDSP). 

The data structure used in our ejection algorithm is a pair (vr, h), where tt is a permutation 
of all the requests R and h G {1,2, ... ,n} is a position in this permutation called hole. In 
order to calculate the objective ^(vr, h) of the solution (tt, h), we choose R* = R, order the 
element of R* according to the permutation tt and then remove the hth element from R*. 
Then, by applying the greedy scheduling algorithm, we obtain a schedule and calculate its 
objective as defined above. 

The basic move in our local search algorithm is swapping of an element with the 'hole', 
i.e., swapping 7r(/i) with and updating h with i for some i ^ h & {1,2, .. . ,n}. There 
are n — 1 options for this move, and we select the one that produces the best solution (vr, h). 
Every time we obtain a new solution {'K,h), we also evaluate the schedule S^^. If S"^ is better 
than the best schedule known so far, we save it. For details of our ejection chain algorithm 
see Algorithms [1] and [21 

Note that, as it was mentioned above, this algorithm does not take care of the urgent 
requests. To make sure that the urgent requests are scheduled as early as possible, we start 
from running the ejection chain algorithm for the urgent requests only (see Section [3]). This 
guarantees that they are scheduled as early as possibl^. Then we fix these requests and run 
the ejection chain algorithm for the regular requests. 

4 Data analysis and generation of test instances 

The algorithms developed in this paper were tested using real data. We used 13 low density 
instances collected between October 25 and December 13 of 2011, each containing approxi- 
mately 100 requests per planning horizon. We also used 31 high density instances collected 

^In fact, we cannot guarantee that the urgent requests are scheduled as early as possible because we do not 
guarantee the optimality of our solutions. However, since the number of urgent requests is usually relatively 
small, corresponding subproblems are easy to solve, and, thus, we believe that the schedule containing only 
urgent requests is close to optimality. 
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Algorithm 1 Ejection Chain Algorithm. 
Input: Requests R. 

Input: Priorities Pr and time window durations Wr for each r & R. 
Input: Objective functions 0(7r) and 4){7r,h). 

Create a permutation tt of i?2, • • • , Rn) such that p7r(«) > Pnii+i) (or p7r(i) = p,r(i+i) and 

[w^ii]] < [wn{i+i)\: or p„(i) = Pn(i+i): [wn(i)\ = L^7r(j+i)J and < for every 

i = 1,2, . . . ,n — 1. 

Save TT to a permutation cr. 

Initialize the idle counter c with n. c keeps track of the number of non- improving iterations. 
Initialize the hole position h with 1. 

While there arc no more then n consecutive non-improving iterations: 
while c > do 

if improvement{Ti , h, a, 10) = 1 then 
Reset the idle counter c with n. 
Save TT to a. 
else 

Reduce the idle counter c by one. 
end if 

Go to the next hole position: if /i = n set ■<— 1 otherwise h -ir- h + 1. 
end while 
return a. 



in August 2011, each containing approximately 300 requests. There are about 10 ground 
stations involved in each instance and these ground stations are typically visible during the 
planning horizon from 4 to 10 times. 

For the low density data, two solutions were given to us: machine- scheduled and human- 
rescheduled. The machine scheduled solutions are generated by the algorithm currently in 
use for mission planning. Such a solution may be infeasible and the human-rescheduled 
solution modifies the machine scheduled solution to remove violations and to satisfy some 
additional considerations known to the operators at that time. Such a data (and solution) 
may be imprecise. For example, our model uses crisp visibility mask boundaries but a 
human operator may use his/her judgement to schedule a job even if a job goes beyond the 
prescribed visibility mask by very small amounts but such a solution will be infeasible as per 
our model. 

The high density instances were given to us as a single file containing data for the month of 
August 2011. In particular, we were provided with a table of all the requests and the schedule 
generated by human intervention. Since we were not provided any explicit instances for each 
planning horizon of 24 hours, we used the following rules to construct the instances. Let Tj 
be the planning horizon of the ith instance, i — 1, 2, . . . , 31. 

1. If a request r is present in the provided schedule, we assign it to the ith instance, where 
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Algorithm 2 Ejection Chain Algorithm: improvement{7r,h,a,d). 
Input: Permutation ir with and a hole position h. 
Input: The best solution a found so far. 

Input: The remaining search depth d, i.e., potentially the longest ejection chain 
if d — then 

return 0. 
end if 

Save the objective value of the partial solution 0(7r, h) to (pp. 
Initialize j with 0. 

for i ■<— 1, 2, . . . , — 1, ^ + 1, . . . , n do 
Swap 7r(/i) with 7r(i). 
if (f){Tr) > (l){a) then 

Save the solution tt to a. 

return 1. 
end if 

if (^(tt, i) > 0p then 

Update j with i. 

Update 0p with 0(7r,i). 
end if 

Swap n{h) with 7r(i). 
end for 
if J 7^ then 

Apply the modification by swapping 7r(/i) and 7r(j). 

if improvement{7r,j,a,d— 1) = 1 then 
return 1. 

end if 
end if 
return 0. 



2. If an urgent request r is withdrawn in the provided schedule, we assign it to the first 
instance where it can be downlinked. 

3. If a regular request r is withdrawn in the provided schedule, we assign it to the zth 
instance, where i is chosen such that the length of the interval Tj fl [6^,6^ — d,.] is 
maximized over i G {1,2,..., 31} (if there are several such is, we chose the smallest 
one). 

Several of the 31 instances were excluded from our test bed because of a significant 
amount of violations found in the provided solution. These are instances corresponding to 
Aug 7, Aug 8 and Aug 31. 
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5 Computational Experiments 



The ejection chain algorithm is coded in C++ and tested on a PC with Intel 17-2600 CPU. 
The low and high density instances as discussed above were used as our test bed. We 
compared our results with the solutions provided by our industrial partner. The summary 
of our experimental results on low density instances is given in Table [H 
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33 


1.1 


1.1 


0.0 











1 


1 





0.3 


1.3 


0.0 


201 


202 


119 


Nov 23 


90 


30 


1.2 


1.1 


0.0 











1 


1 





0.4 


3.1 


0.0 


371 


379 


79 


Nov 29 


108 


31 


1.0 


1.0 


0.0 











1 


1 





4.5 


4.5 


0.0 


234 


234 


102 


Dec 1 


110 


28 


1.9 


1.6 


0.0 











2 


1 





5.7 


11.3 


0.0 


536 


975 


237 


Dec 2 


134 


27 


0.3 


0.3 


0.0 




















0.4 


0.4 


0.0 


564 


564 


65 


Dec 7 


108 


36 


1.1 


1.1 


0.0 











1 


1 





1.4 


0.5 


0.0 


573 


575 


151 


Dec 13 


94 


24 


2.1 


2.1 


0.0 











2 


2 





0.2 


0.2 


0.0 


566 


566 


40 


Average 


107.6 


30.0 


1.5 


1.5 


0.0 


0.0 


0.0 


0.0 


1.3 


1.1 


0.0 


1.8 


2.0 


0.0 


617 


704 


162 



Table 1: Comparison of the given machine- scheduled solutions (M-S), human- rescheduled 
solutions (H-R) and the results of our ejection chain algorithm (EC) for the low density 
instances. 

Each row of the table corresponds to one problem instance. For each of the algorithms 
and for each of the instances we provided how much worse is the solution comparing to 
the best known solution for this instance ('To the best, %'), the number of unscheduled 
urgent requests ('Urg. unsch.'), the total number of unscheduled requests ('Total unsch.'), 
the average delay of the urgent downlinks ('Avg. urg. del., s') and the average delay of the 
regular downlinks ('Avg. reg. del., s'). Since the machine-scheduled solutions are sometimes 
infeasible, we consider a downlink to be scheduled if it falls into the planning horizon and is 
not present in the violations table provided with the machine- scheduled solution. The delay 
of a request r is measured where is the earliest time when r can be downlinked. 

Note that the comparison of our schedules to the human-rescheduled solutions may be 
incorrect. Indeed, we may not know some of the constraints that are taken into consideration 
by the operator. Thus, we focus on the comparison to the machine-scheduled solutions. 

The computational results show that the solutions produced by our algorithm dominate 
the provided machine-scheduled solutions for all of the instances. In particular, our algorithm 
always schedules all the downlinks while usually there were several unscheduled requests in 
the given solutions. Our algorithm never delays the urgent requests while it sometimes 
happens in the given machine-scheduled solutions, i.e., there are large delays (around 3 
minutes) for the Nov 17 and Nov 29 instances. In each case, the delay is mainly caused by 
only one downlink that the algorithm failed to schedule to the beginning of the time window. 

Note that, since our algorithm schedules all the requests for every low density instance. 
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the optimality of schedules does not depend on the value of the constant a (see Section |2]). 

The running time of our algorithm is always very reasonable and never exceeds 10 seconds 
in our experiments. 

Computational results for high density instances are given in Table |2l 



Instance 


n 


Urg. 


To the best, % 


Urg. 


unsch. 


Total unsch. 


Avg. urg. 


del., s 


Avg. rcg 


. del., s 


H-R 


EC 


H-R 


EC 


H-R 


EC 


H-R 


EC 


H-R 


EC 


Aug 1 


211 


34 


13.5 


0.0 








53 


19 


0.3 


0.0 


2013 


457 


Aug 2 


301 


35 


20.5 


0.0 








123 


67 


0.6 


0.0 


3788 


2725 


Aug 3 


324 


42 


20.7 


0.0 








120 


48 


0.5 


0.0 


2329 


2248 


Aug 4 


294 


41 


21.1 


0.0 


2 





116 


52 


0.5 


0.0 


2025 


1248 


Aug 5 


260 


36 


13.3 


0.0 








85 


49 


0.5 


0.0 


2621 


1021 


Aug 6 


356 


46 


20.0 


0.0 


2 


1 


130 


61 


0.4 


0.0 


2555 


1563 


Aug 9 


259 


38 


18.7 


0.0 


3 





87 


30 


92.4 


0.0 


2448 


1867 


Aug 10 


304 


54 


20.9 


0.0 


1 





98 


33 


0.4 


0.0 


2230 


1338 


Aug 11 


278 


39 


19.4 


0.0 


3 





90 


36 


0.5 


0.5 


2463 


904 


Aug 12 


230 


25 


16.8 


0.0 








80 


28 


0.5 


0.0 


2880 


2096 


Aug 13 


324 


41 


16.7 


0.0 


1 





101 


38 


0.4 


0.0 


2045 


868 


Aug 14 


285 


34 


18.8 


0.0 


4 


3 


104 


54 


0.4 


0.0 


1558 


818 


Aug 15 


259 


30 


18.6 


0.0 








93 


43 


0.4 


0.0 


1692 


1303 


Aug 16 


299 


39 


29.5 


0.0 


4 





139 


54 


0.4 


0.0 


1374 


1849 


Aug 17 


303 


30 


23.0 


0.0 








128 


65 


0.5 


0.0 


2664 


1193 


Aug 18 


275 


33 


20.6 


0.0 


2 


1 


115 


54 


0.4 


0.0 


2557 


1819 


Aug 19 


254 


30 


18.8 


0.0 








93 


34 


0.4 


0.0 


2210 


1817 


Aug 20 


295 


39 


19.4 


0.0 








96 


35 


0.4 


0.0 


2768 


1124 


Aug 21 


272 


24 


18.7 


0.0 








100 


41 


0.4 


0.0 


2531 


1479 


Aug 22 


253 


36 


18.0 


0.0 








91 


45 


0.5 


0.0 


1737 


1486 


Aug 23 


307 


44 


18.6 


0.0 








111 


45 


0.3 


0.0 


2450 


1421 


Aug 24 


300 


40 


22.9 


0.0 


1 





122 


56 


0.4 


0.0 


1465 


1146 


Aug 25 


269 


40 


36.2 


0.0 


11 





112 


48 


2208.0 


0.0 


6051 


1664 


Aug 26 


301 


37 


19.2 


0.0 








116 


44 


0.5 


0.0 


1955 


2152 


Aug 27 


282 


40 


18.5 


0.0 








89 


33 


0.4 


0.0 


2720 


921 


Aug 28 


282 


34 


21.9 


0.0 


5 





108 


41 


0.5 


8.0 


2765 


1014 


Aug 29 


267 


36 


20.5 


0.0 


3 


1 


95 


46 


0.4 


0.0 


2991 


1393 


Aug 30 


315 


44 


21.0 


0.0 


1 





116 


46 


0.3 


0.0 


2081 


1582 


Average 


284.3 


37.2 


20.2 


0.0 


1.5 


0.2 


104.0 


44.5 


82.6 


0.3 


2463 


1447 



Table 2: Comparison of the human- rescheduled (H-R) solutions with the results of our 
ejection chain algorithm (EC) for the high density instances. 

Unlike it was for the low density instances, our ejection chain algorithm is not able to 
schedule all the requests in the high density instances. However, it schedules 60 requests 
more, on average, than in the given machine-scheduled solutions. It also schedules more 
urgent requests and provides a smaller average downlink delay. The running time of the 
algorithm is below 4 minutes on average and it never exceeds 10 minutes. 

Another experiment that we conducted was to schedule as many requests from the high 
density instances as possible. For this purpose we defined the instances as follows. The 
Aug 1 instance includes all the requests that can be downlinked on August 1. The Aug 2 
instance includes all the requests that can be downlinked on August 2 except the requests 
scheduled for Aug 1. The Aug 3 instance includes all the requests that can be downlinked 
on August 3 except the requests scheduled for Aug 1 and Aug 2, etc. 
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We call such instances rolling. Note that rolling instances may be different for different 
algorithms and, thus, comparison of the performance of the algorithms for every particular 
instance is meaningless. 

The results of our experiment with the rolling instances are reported in Table [31 



Algorithm 


Objective 


Scheduled 


Unscheduled 


Avg. delay, sec 


Running time, hours 


Human- rescheduled 


219.1 


5056 


3428 


9486.6 




Construction 


293.8 


6513 


1971 


12286.6 


0.0 


Ejection chain 


304.1 


7084 


1400 


8995.6 


3.9 



Table 3: Comparison of the given human- rescheduled solution with the results of our con- 
struction heuristic and our ejection chain algorithm for the rolling instances. 



It takes around 4 hours for our ejection chain algorithm to schedule 7086 out of 8484 
requests for the entire month of August, compared to 5056 requests scheduled in the human- 
rescheduled solution. Note that our algorithm also improved the average delay and, hence, 
the quality of the solution. 

Our construction heuristic proceeds in almost no time and schedules more requests than 
in the given human-rescheduled solutions. However, the average delay in this solution is 
higher. 

Although such an experiment may not reflect the real procedure of scheduling downlink 
requests, it shows the potential for improving downlink schedules. 

5.1 The parameter a 

Recall that the objective function for the SIDSP uses a parameter a that represents the 
relative preference of the downlink delays over unscheduled requests. Our ejection chain 
algorithm was able to schedule all the requests for the low density instances and, hence, the 
solutions obtained for these instances were almost independent of the value of a. In contrast, 
for the high density instances, careful selection of the value of a was crucial in achieving high 
quality solutions. In Figure [1] we show how the performance of our ejection chain algorithm 
depends on the value of a for the high density instances. 

The tendency is that as a increases, the delay in scheduling a request decreases and the 
number of unscheduled requests increases. Hence, a can be used as a parameter to adjust 
the trade-off between the number of scheduled requests and the average delay. 

Note that the performance of the algorithm is poor for a = 0. This is because the objec- 
tive in this case does not depend on the delays making the search landscape more discrete 
and, thus, complicated for optimization. It is also worth mentioning that the difference in 
the solutions obtained for different values of a is relatively small. 
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Figure 1: How the value of a influences the performance of our ejection chain algorithm. 
The experiments were conducted for the high density instances. 



6 Conclusions 



In this project, we formalized the satellite downlink scheduling problem and proposed an 
efficient heuristic solution approach. Computational experiments show that our algorithm 
clearly outperforms the currently implemented algorithm, and the running time of our algo- 
rithm is practically negligible to reasonable. 

It is worth noting that we understand that the proposed problem definition may not be 
precise and some additional constraints may be found out later. Our approach, however, is 
very flexible and we believe that it can easily incorporate potential additional constraints. 
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